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_ . Abstract. We study the limit of asymptotically free massive integrable models in 

C J ' which the algebra of nonlocal charges turns into affine algebra. The form factors of 

fields in that limit are described by KZ equations on level 0. We show the limit to 
be connected with finite-gap integration of classical integrable equations. 
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fN^ . 1. Motivations. 

l/~j I The nonlocal symmetries of integrable models of quantum field theory in two 

^D i dimensions were first studied several years ago [1]. The reason for that being in 

^— ^ I attempts to understand the quantization of asymptotically free models. Being al- 

most forgotten for some time the nonlocal symmetries returned to the field rather 
Qs I indirectly, namely, through the finite-dimensional quantum group symmetries of 

CFT. Now it is understood that the integrability of massive models is closely con- 
nected with possessing infinite-dimensional algebra of nonlocal symmetries [2-6]. 
The local integrals constitute a center of it. The algebra of nonlocal symmetries is 
qj I always a Hopf algebra, particles transform under its finite-dimensional representa- 

<^ . tions while quasilocal fields constitute infinite-dimensional multiplets with highest 

vectors corresponding to the local fields. The S-matrix is nothing but universal 
R-matrix specified onto finite-dimensional representations while the braiding of the 
K% ' multiplets of quasilocal fields is described by the universal R-matrix specified onto 

^ . the tensor product of two Verma modules. Moreover, it was shown [6] that the form 

factors which put together particles and fields can be considered as solutions of de- 
formed Knizhnik-Zamolodchikov equations [7] which is quite natural since these 
equations use to relate finite-dimensional representations of deformed loop algebras 
with infinite-dimensional highest weight representations. 

A question which should be asked is the following. Suppose we have an inte- 
grable model with symmetry under certain deformation of the loop algebra. Then 
what is the meaning of the "classical limit" which moves the deformed algebra into 
undeformed one. The answer to this question is not trivial as we shall see. For 
many models this limit does not look to make much sense. The point is that the 
limit is related not to the rescaling of space coordinates (or momenta) as conformal 
limit is, but rather to rescaling of rapidities of particles which are logarithms of 
momenta. In a sense the limit in question is even opposite to conformal limit in 
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which we throw away all the logarithms preserving powers, here we preserve loga- 
rithms, but throw away powers considering the asymptotic expansions. However, 
the limit seems to be reasonable for asymptotically free models. So, we return to 
the same idea as years ago: to use nonlocal symmetries in asymptotically free the- 
ories, certainly we return to this idea with new experience. It should be said also 
that there exist certain problems with correspondence between local integrals of 
motion for classical and quantum cases in the theory of asymptotically free models. 
The reason for that is very fundamental, that is why our consideration will be not 
Hamiltonian one. 

Particular (and, probably, personal) interest in this business is the following. We 
have very nice formulae for final physical quantities in both classical and quantum 
theory of integrable models. However existing way of quantization [8] (which has 
clear Hamiltonian meaning) does not allow to proceed directly from ones to others, 
it starts with beautiful things (R-matrix, algebraic Bethe anzatz), but then leads 
through usual jungle with rcnormalizations and all that. We would like to have 
direct way even if it has less clear meaning. Similar ideas were applied to study of 
correlation functions of vertex models of statistical mechanics in recent works from 
Kyoto school [9]. 

Let us consider a particular asymptotically free model. Namely, let it be the 
su{2) chiral Gross-Neveu model. The spectrum of the model contains one two- 
component particle. The two-particle S-matrix is given by [10]: 

V ZTTl I \ 2. 2-Kl / 

where Pi^2 is permutation. If we are interested in local operators which transform 
under integer-spin representation of global SL{2) isotopic group all of them can be 
obtained as descendents with respect to different conservation laws of one nonlocal 
operator called ^{xo,xi). The connection with important local operators will be 
explained later. The operator 2t(xo,a:i) is scalar with respect to both isotopic and 
Lorenz transformations. In form factor approach [11] it can be presented as follows: 



oo _. 

2l(xo,a:;i) = T^ — r- 



mini 

m,n— 



dai ■ ■■ dujn / dPi-- ■ d(in Z*^ (ai) • • • Z*^ (a™) 

) -'— oo "'— oo "'— oo 

X F(a„, • • • ,ai|/3i, • • • , /^n)^":::;^': ^'"(/3n) • --Z'^Pr) 

X exp(^ip^(a,)a;^ -^ip^{l3i)Xf,) (1.1) 

where Z* , Z are Zamolodchikov-Faddeev creation-annihilation operators of parti- 
cles, p^{a) is one-particle energy-momentum (p^{a) — M(exp(a) + (— l)''exp(a))), 
e, (5 = ± are isotopic indices. F are form factors. They satisfy the crossing symme- 
try requirement: 

F(a„,-- • ,q;i|/3i,-- • ,PnTi^^^'\t^ = cs^,s[ ■■ •c^„,<5; 
X F{am,- ■■ ,ai,/3i +7r2 -iO, •• • ,/3„ + ttz - iO)'"'"' '^i'"^!'"' ^''", c = ia"^ , 
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It should be mentioned that for the operator 21 the two particle form factor has a 
double pole at /32 = /3i + ttj the corresponding singular part being normalized as 
i(32-f3i-Tri)-^c'^^'-\ 

The model allows abelian and non-abelian symmetries. The first are the local 
conservation laws with all possible odd spins. The one particle eigenvalues of the 
local integrals Is {s — 2m+ 1, m — — oo, • • • , oo) are equal to exp(sQ;) the spectrum 
of them being additive which can be expressed in more formal way by the following 
comultiplication formula: 

A(/s) = Is®l + l®l8 

The nonlocal symmetries together with S'i(2)-charges constitute the Yangian {Y) 
[2]. The generators of Y are J^jJf, a = 1,2,3 whose action onto one particle is 
described by cr", acr" respectively and the comultiplication is given by 

A(Jo") = Jo" ® 1 + 1 ® Jo", 

A(jn = Jf (» 1 + 1 ® Jf + 7:ir^^4 ® Jo" (1.2) 

Every operator acting in the space of states (Hp) which is the Fock space of par- 
ticles has descendents with respect to both local and non-local symmetries defined 
through ajoint action. One can think of the space Hp as of 

oc „ 

^p=0 J ^ft®---®^&„ (1-3) 

"=l/3i<---</32„ 

where Vp- is the space of the representation p/^.. The operators Z*{/3), Z{/3) ei- 
ther add or remove one space Vg. In particular the energy- momentum tensor and 
currents are the following descendents of the operator 21: 

Tf_,,i^{xo,xi) = efj^^^>e^yPf^'P^>{^{xo,xi)), j°(xo,a;i) = e^,^'P^' Ji"(2t(xo, xi)) 

where Pq = Ii+ I^i, Pi = /i — /_i. There are more general relations. For example 
the density of the local conservation law /„ {ds(xo, xi)) is obtained in the following 
way: 

ds{xo,xi) = PoIs{^{xo,xi)) 

The last relation is very much in common with the classical relation identifying the 
densities of local conservation laws with the derivatives of logarithm of r-function. 
So, it is reasonable to think of the operator 2l(a;o,a:i) as of the quantum analog of 
log(T). Certainly, this is a formal analogy and in order to make it more instructive 
we have to develop the relation with the classical r-function. This is the main goal 
of the present paper. For the reasons which will be explained later the connection 
in question can hardly be achieved in Hamiltonian formalism, so, our strategy will 
be to find a "good" formula in classics in which all the objects involved in quan- 
tum formula (the operator 21, the form factors, Zamolodchikov-Faddeev operators) 
will have their classical analogs. This is what we are going to do in that paper. 
But before proceeding in that direction we have to remind certain facts about the 
symmetries of the quantum model. 
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2. The Yangian symmetry. 

In this section the brief account of the recent developments concerning the dy- 
namical symmetries of integrable models will be presented. 

As it has been already mentioned the model under consideration allows nonlocal 
conservation laws Jf which together with isotopic charges Jq constitute the infi- 
nite dimensional algebra called Yangian [12]. It is explained in [5,6] that for the 
application to the description of local and quasilocal operators in the theory one 
has to add one more generator {Jti) to the algebra. The full algebra generated by 
Jq , J°, J'^i is called the double of Yangian being denoted by D{Y). One can think 

of D{Y) as of a deformation of the afline algebra sl{2) with natural identification 
of the generators. This will be explained more precisely later. The Yangian y is a 
subalgebra of D{Y) which acting onto the fields in the theory creates descendents. 
The introducing of the second half of the algebra was originally motivated by the 
problem of description of the commutation relations among these descendents [5]. 
The commutation relations are described as follows. Consider a local operator 
(p{xQ, xi) which can be for the simplicity taken as invariant under the isotopic alge- 
bra. Now create all possible descendents of the operator acting on it by arbitrary 
number of the operators Jf: 

^{xo, xir^- ■"" = J^' ■ ■ ■ Jf"((^(xo, a;i)) 

We can combine these operators into <i> considering ai, • • • , a„ as multiindex. So, 
$ is an infinite column of operators which can be considered as belonging to 
End{Hp) (g) W where Hp is the space of states, W is the Verma module created by 
action of the operator J". Now, if we consider two towers of operators ^(xq,xi) 
and ^{xq,x'i) all the products of their elements can be combined in the following 
object: ^i{xo,xi)^2{xa,x[) where two copies of W {Wi,W2) are considered and 
$1 G End{Hp) (g) Wi, $2 € End{Hp) (g) W2- The commutation relation in question 
can be written down in the following way [5] : 

$i(xo,a;i)$2(a;o,a;i) = Ri^2^i(xt^, x^)^2{xn, xi), for xi > x[ 

Here Ri_2 is D(Y) R-matrix acting in the tensor product of two Verma moduli. 
That means that in these relations we effectively consider the local field (p as that 
annihilated by the operator J"]^. 

It is explained in the papers [6] that the situation can be inversed. We consid- 
ered the space of particles and came to the conclusion that the local and quasilocal 
operators acting in the space behave, as far as the commutation relations are con- 
sidered, as being combined into Verma moduli. The local operators themselves 
are identified with the highest vectors of the Verma moduli. We can consider now 
the space of fields instead of the space of particles. The space of fields (Hf) is a 
suitable collection of Verma moduli. Particular local field (pi is identified with a 
highest vector |0)i satisfying the requirement 

J-i|0). -0 

We consider also a dual vacuum (0| which is defined by the relations 

(0| Jo" - 0, {0\J^ = 
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The algebra D{A) allows beside of the infinite dimensional representations the 
finite-dimensional ones. We will be interested in two-dimensional representation 
(pis) depending on a parameter /3, the detailed description is given in [6]. One can 
define the vertex operators V^IP) which belongs to End{Hf) (g) C^ (e = ± is C^ 
index). Acting in Hf the operator V'^{P) transforms under the ajoint representation 
with respect to the representation pfj: 

Then the vacuum expectations of the operators V^ (/3) given by 

satisfy the deformed Knizhnik-Zamolodchikov equations, in usual tensor notation 
they take a form: 

(0|F(/3i) • • • V{pj + 27r) • • • V{p2nm^ = 

= S2n,jiP2n - Pj - StTi) • • • Sj + i^j{Pj+i - f3j - 27ri) 

x(om/3i)---n/?j)---n/52„)io), (2.1) 

As it has been explained in [12] these equations are consistent with the symmetry 
property 

5m+i(A-A+i)(0|F(/3i),--- ,F(A),T^(/3.+i),--- ,n/?2„)|0). - 
= P,,,+i{0\V{pi),- ■ ■ ,V{p,+i),Vm,- ■ ■ ,V{(32nm^ (2.2) 

Being supplied with the symmetry property (2.2) the equations (2.1) appear to be 
the same as the basic requirements for the form factors. There is an additional 
equation on residues in the form factor bootstrap approach [11] which we do not 
present here. This equation can be interpreted as a form of operator product 
requirement [6]. It is responsible for the special choice of blocks of the vertex 
operators. Namely, the following remarkable identification can be done 

FiPu- ■ ■ , I32nr'- '''- ^ {0\V'' (/3i) • • • T/^- (/32„)|0) (2.3) 

where |0) is the highest vector of spin zero. The consequence of the Verma moduli 
in the RHS is taken as follows (we indicate the spins of the highest vectors) 

0^1/2^1^1/2^1^ > 1/2^0 (2.4) 

We put the asymptotical equivalence in (2.3) for the reasons explained in [6], we 
shall briefly explain the point later. The form factor can be considered as a matrix 
element of the operator 21 in the space of particles. More generally the relation 
(2.3) can be rewritten as follows 

(mc|Z^n/3i) • • • Z'--{|32n)^^m^ac) ^ (0|y^H/3i) ' • • "l^''"(/32„)|0). 
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We do not expect other local operators then the descendents under the action of 
local integrals on 21, and on the currents j" to exist in the theory. There is also 
paraferniionic field of spin j with respect to Lorentz transformations (kink field) . 
From the point of view of the representation theory these three fields correspond 
to spin 0,1,1/2 highest vectors respectively. The sequence of the Verma moduli in 
form factors of these operators is the same as in (2.4) the only difference being 
that the sequence for currents (kinks) terminates at spin 1 (1/2) representation. 
This is a variant of "rationality" of the theory. The reason for the phenomenon to 
occur supposed to be in the absence of mixing of different solution of deformed KZ 
equation by braiding which is discussed in details in the papers [6]. 

Now we are going to explain the asymptotical character of the equivalence in 
the relations (2.3). This matter was considered in [6,12]. The situation is briefly 
as follows. In LHS of (2.3) we have the analytical solutions of (2.1) while the 
RHS contains the vacuum expectation of £'(y)-vertex operators. The RHS can 
be in principle calculated directly (without use of KZ) using the definition of the 
vertex operators. It is quite clear that the only thing we are able to get in that 
way are certain power series in /3's. But as it follows from the explicit formulae 
for the analytical solutions they are transcendental functions of /3's (the explicit 
formulae are given in [11]); roughly, they have the properties of F-functions. The 
only reasonable connection between these two types of objects is asymptotical one. 
One can have in mind the following analogy: F-function (RHS) and its asymptotical 
series (LHS) both satisfy the functional equation for the F-function (dKZ-equation). 
The difference between the function and its asymptotics in the case is due to the 
exponential in (3 contributions which are considered as the contribution due to the 
intermediate state created by local integrals which should be added to the D{Y) 
[6]. 



3. The classical limit. 

As it has been already said the Yangian double is a deformation of the algebra 
sl{2). Let us ignore for the moment the physical content, and consider the formal 
aspects of the limit. The typical relation for us will be (1.2). By rescaling of Jf 
the relation can be rewritten as 

A(Jo ) - Jo ® 1 + 1 «) Jo , 

A(Jf ) = J^* ® 1 + 1 ® Ji° + -^/"^"^ <» -^0 (3-1) 

where h is arbitrary constant. Certainly , the generator J^j^ should be rescaled as 
well. Now consider h as Planck constant and take the limit h -^ +iO ( we prefer 
considering imaginary h, so, probably, temperature is better analog for it then 
Planck constant). In the limit the equations (3.1) turn into trivial comultiplication 
formulae for the generators of sl{2), all other I?(y)-relations turn into s^(2)-ones 
in the limit. Also the dKZ equations turn formally into usual sl(2) KZ equations 
on level zero. It is clear that the rescaling of the generators we did for the dKZ 
equations is equivalent to the rescaling of the rapidities: Pi = ^A^ , h -^ +iO 
while Xi arc fixed. So, formally for the functions 

f/\ \ ^ r^fh\T?f'^'^'^'^^ 27riA2„, 



{C{h) is a normalization constant) the equations (2.1) turn in this limit into 

'A. 
'dX, 



(jT+T. ^^^^ (^* - ^j))/(Ai, • • • , A2„) = 



where r is the classical r-matrix: 

rij(A, - Aj)) = 



erf (g) a"} 
Xi — Xj 



The connection between deformed and undcformcd KZ equation is described ex- 
plicitly in [13]. Let us outline the basic points. 

First, the dKZ equation allows the same number of solutions as undeformed 
one. The solution corresponding to the form factors is a special one. Its particular 
character is explained by (2.4), we shall also describe it from other point of view 
soon. 

Second, the solutions are in one-to-one correspondence in the following sense. 
Consider some solution of dKZ for /3i < /32 < • • • < f32n then its asymptotics 
/3j = ^Ai , ft- ^ +iO is described by a solution of KZ. In the paper [13] the 
explicit formulae are presented for solution of dKZ which correspond to all possible 
solutions of KZ in that sense. 

Third, there is essential difference in the properties of the solution of dKZ in 
comparison with those of KZ: braiding does not mix different solutions of dKZ. 
The difference is explained by the asymptotical character of the correspondence 
above. Braiding does not commute with taking the asymptotics. 

The third point mentioned makes a real difference between the deformed and 
undcrformed case. In particular, in the deformed case it makes sense to consider 
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the form factors of 21 themselves without consideration other solutions to dKZ 
equations. But taking the asymptotics of the the form factors we have to consider 
not only the solution of KZ we obtain, but also those which are connected with it 
by braiding. We shall call this part of solutions of KZ the main part. 

We do not give here the explicit formulae for the solutions of dKZ which can 
be found in [11,13], but we have to write down explicitly the solutions of KZ. 
These formulae concerns the particular case of zero central extension which causes 
essential simplification with respect to the general case. So, let us explain first how 
the simplifications appear. In Varchenko-Shcchtman [14] formulae for the solutions 
of KZ one has the following structure: 

Y[iX, - Aj-)5(^ / dn • • • dTk J|(t„ - a,) <^ n (^' " T-m)^^ 
i<j m,i l<7n 

xi?(Ti,--- ;Ai,---) (3.2) 

where r's are integration variables (corresponding to screenings), _R is a rational 
function of its variables. For k = the exponents in (3.2) turn into 1/4, 1/2, 1 
respectively. The first one is not so essential, the second causes the hypercllip- 
tic character of the integrals, the third causes the absence of branching between 
the integration variables which simplifies essentially the choice of contours of in- 
tegration. Also one can present n;<m('''' ~ ''m) ^^ Vandcrmond determinant, and 
rewriting properly the function R perform the integration over the columns of the 
determinant getting a determinant of single integrals instead of the niulti- integral 
of determinant. Let us present the result to which these manipulations should lead 
(which actually has not been obtained in that way, but directly from the classical 
limit the solutions of dKZ solutions). 

Consider the components of /(Ai,--- ,A2„): /(Ai,--- ,\2nY^'"'''^^"- We are 
looking for the singlet solutions of the equations which means in particular that 
^ei = 0. For each particular component /(Ai,--- , X2nY^' ''^^" the multiindex 
ei, • • • , ^2n induces a partition of _B = {1, • • • , 2n} into T ~ {ik}k=i • ^it — + ^^'^ 
T' = B\T = {jk}k=i ■ ^jk = ^- Different solutions will be parametrized by the 
sets 7i, • • • ,7n-i which will be specified later. The solutions look as follows 

/^,,..,^„.,(Ai,...,A2„)'--'^-=n(^«-^j)" 

i<j 

X n (A.-A,)-Met|| / 0(T|T|T')dT||,x3 (3.3) 

ieT.jeT' "^T' 

C,j are the following differentials on the hyper-elliptic surface (HES) w^ ~ P{t) = 
n(T-A,): 

S ar T ^, ar T ^^^^^ 

where [ ]o means that only the polynomial part of the expression in brackets is 
taken. The differentials Q are of the second kind: they have singularities at oo*, 
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but their residues at the infinities are equal to zero. It should be mentioned that 
the singular part of the differential Q is independent of the partition of A: 

Cj(T|r|r') - Q{t\Ti\T[) = of the first kind 



First kind differentials in our case are of the type: aj ~ t^^^ j ^P(t), \ < i < n—1. 
The contours 71, • • • ,7n-i are arbitrary cycles on the HES ( notice that its genus 
g equals n — 1 ). 

It can be shown that the asymptotics of the form factor F corresponds to the fol- 
lowing special choice of 71 , • • • , 7n- 1 : they are taken as canonical a-cycles ai , • • • ,ag 
where the cycle a^ surrounds the cut between A2i_i and A2i, i — 1, • • • ,n — 1 (we 
suppose A's to be ordered: Ai < A2 < • • • < A2„). 

Now, let us turn to the problem of braiding. The equations (5.3) are invariant 
under the permutation Aj <-^ A^ and simultaneous permutation of the associated 
spaces. Let us denote the operation of the analytical continuation Ai ^^ A^+i and 
permutation of corresponding spaces by i?i,i+i. Then we are supposed to have a 
formula of the type: 

-Di,i+l/7i,--- ,7„-i (Al, • • • , A2ri) = 

where C are some constants. There is one interpretation of the braiding following 
from the formula (3.3) : Ai <-> Ai+i corresponds to certain element of the modu- 
lar group Sp{2g, Z) of the HES. Clearly under the braiding the determinants (3.3) 
transform under g-th exterior power of the vector representation of the group. It 
can be shown that this interpretation implies that the solution / corresponding to 
7i I ■ ■ ■ 1 7g = 0,1, ■■ ■ T<ig raixes by braiding only with those solutions which corre- 
spond to other choices of half-basis of homologies {g contours on the surface with 
zero intersection numbers). This is the main part of solutions of KZ defined above. 
There is the same number of elements in the main part as the number of indepen- 
dent half-bases, i. e. C|J^2 ^ '^2q+2 of them. On the other hand it is well known 
nowadays that the braiding is described by the finite-dimensional quantum group 
SL{2)q [15], in our case q = —1. The number of solution coincides with the the 
multiplicity of the one-dimensional representation in the tensor product of 2n two- 
dimensional representations ot SL{2)q which is in generic situation CJ„ — C^,7 
(the same as above since g = n — 1). So, it is natural to suppose that the main 
part of solutions for the case q = —1 allows continuation to other q while the rest 
of solutions is special for g = — 1. In what follows we shall need only the main part. 
There is a beautiful way of rewriting the formulae (3.3) for the case of the main 
part of solutions in terms of ^-functions. We shall consider that in the next section. 
Now let us turn to the problem of physical interpretation of the classical limit 
in question. We have to understand what does the limit D{Y) -^ sl{2) mean. In 
that limit, in particular, the antipode-square automorphism of D{Y) turns into 
differentiation for sl{2): 

s^^l + hD (3.5) 

We do not denote D by familiar notation (i-i) to avoid a confusion: D is not 
a derivative in the space-time. There is a strange thing about this limit. The 
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antipode-square (s^) is identified with the rotation of the space-time in the theory 
[6], for two-dimensional representations (particles) it corresponds to the shift of 
rapidity by 2iri. So, the operation, s^ corresponds to, in the space-time of QFT 
model is essentially finite being of topological character, hence the consideration of 
its infinitesimal limit (3.5) does not look very reasonable. However we do will to 
consider the limit because from mathematical point of view the limits of certain 
objects (e.g. form factors) do make sense. Thinking more of this situation one 
comes upon an idea that in the limit in question the very notion of the space-time 
of the QFT model should be lost. Form factors are the objects referring to one point 
in the space-time, so their limits in this strange situation might make sense because 
the notion of one point might remain (it simply appears to be isolated from the rest 
of the world). Certainly form factors also do not survive in full meaning. First, we 
deal not with a real limit, but with asymptotical equivalence. Second, the difference 
disappears between two-dimensional representation and dual to it which contained 
for D{Y) essential crossing shift, and the two-dimensional representation can not 
be considered as those corresponding to particles. Thus, the only reasonable guess 
we can make is that the QFT in the limit splits into a family of systems with finite 
degrees of freedom, the limits of form factors being somehow connected with these 
systems. The situation is difficult also because wc do not expect any connection 
between these finite-dimensional system and QFT on Hamiltonian level. 

Let us consider the situation in opposite direction. Suppose we have certain 
family of classical systems every of which has a finite number of degrees of freedom. 
These systems should be unified through the fact of possessing the symmetry under 
sl{2). Now we perform a quantization of these systems which essentially leads 
only to the quantization of sl{2) (sl{2) -^ D{Y)). It should be mentioned that 
there is a jump in quantization of sl{2): as far as we got a small h in D{Y) it 
can be immediately rescaled to finite one. So, the quantization provides a finite 
operation (s^). Now we introduce a space-time, and identify this operation with 
the rotation of this space-time. The coordinates are introduced as those respecting 
this interpretation. As soon as it is done we get the notion of particle, and one 
space-time point we started with appears to be able to interact with others through 
exchange of particles. This interpretation seems to be reasonable. In particular it 
solves the problem of disagreement between local conservation laws in classical and 
quantum integrable asymptotically free models. In a context close to the present 
consideration this problem was discussed in [6]. 

So, now we have to explain what classical systems with finite number of degrees 
of freedom and s/(2)-invariance we have in mind. The answer is implied by the 
structure of formulae of the solutions of KZ on level zero. They are connected with 
HES. So, it is natural to suppose that the classical systems in question are stationary 
finite-gap solutions [16-17] of classical soliton equation with sl{2) symmetry. These 
systems will be described in Section 5. But before doing that we have to explain 
the connection between the solutions of KZ on level zero and Riemann ^-functions. 
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4.KZ equations on level zero and Riemann theta- functions. 

In this section we shall derive a formula which express the solutions of KZ equa- 
tions on level zero in terms of Riemann ^-functions. Let us first fix the notations 
and introduce necessary definitions. 

We consider the HES S of genus g with 2g + 2 real branching points ordered as 
follows: Ai < • • • < X2g+2- Let us put the cuts between the points X2i-i, X2i for 
2 = 1, ••• ,5 + 1. The a-cycles on the surface are taken in canonical way: the cycle 
ai surrounds the cut X2i-i, \2i, for z = 1, • • • , (/. The &-cycle bi {i ~ 1,- ■ ■ ,g) starts 
from one bank of the cut X2i-i,X2i, reaches the cut A2g+i,A2g+2 by one sheet, 
and then returns th the other bank of A2i-i,A2i by another sheet. There are g 



nonsingular differentials on the HES: ctj — t^ / ^JP{t). The normalized first kind 
differentials u)j are linear combinations of ct's satisfying the condition: 

LUj = 5ij. 

The martix of periods 17 is defined as 

bi 

The period matrix is symmetric due to Riemann bilinear identity. 

The second kind differentials possess singularities, but their residues at the sin- 
gular points vanish. There is one type of the second kind differentials of particular 
importance. The differentials of this type are obtained by erasing the dependence 
on one argument of the two-differentials (^^{x, y) defined on S x E. The differential 
uP'{x,y) possesses the only singularity at the diagonal being normalized as 

w^(a;, y) = r:^dxdy, a; - y 

[x - yY 

Being considered as a differential in one variable (say x) it satisfy the normalization 
condition: 

^2 ^ 0, Vi 



Finally, the differential oj"^ is symmetric: bj'^(x,y) — uY(y,x). 
Riemann ^-function is defined as follows: 



j|r2) = 2, cxp{7ri ?7i*r2TO + 27rzz*m} 



mGZs 

where z € C^ . The periodicity property says 

e{z + A' + fiA"|0) = exp{-7ra" *0A" - 27riA"z}6'(z|0), A', A" e Z^ 

The ^-function with characteristics is defined by: 

e[r]]{z\Vt) = exp{7ri77" *rjr/" + 27ri(z + r]'Y'q"]e{z + r/' + VLri"\n) 

where rj — {rj' ,ri") is a characteristic: the vectors 77', r?" belong to R^ [R is the field 
of rational numbers). We will be interested in the case when 77 is a half- period 
characteristics which means that rj' ^-q" G 1/2Z^. 
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The definition of ^-function refers to the particular choice of a and b cycles. 
However, the ^-functions defined with respect to different choices of the homology 
bases are connected due to the modular property of ^-function. Suppose that we 
have two homology bases connected via a transformation from ri,2 G Sp{2g, Z): 

B')-it :){b)- it :)-^p^^-^>-Zi-o<^^^ 

"diag" means the vector composed of the diagonal entries of the matrix. The 
transformation law for ^-functions says that [19]: 

0K](z'|f7') =7(dct(Af))l/2exp{^^^z,z,^i^^^}0M(z|f7) (4.1) 

where M ^ {eft + d), z = M*z', f7' = (aO + 6)(cO + rf)"\ the characteristics 
changes as follows 

C) ^\-b d ) \ij" 

7 is irrelevant for our goals 8-th root of unity. 

Let us denote the set {1, 2, • • • ,2(7 + 2} counting the branching points of S by B. 
The most significant property of HES is that they allow one-to-one correspondence 
between half-periods and the subsets of B with even number of elements (mod 
identification T ~ B\T) [19]. The correspondence is achieved as follows. Consider 
the subset T (#r = 2m) and divide it into two subsets Ti,T2 such that #Ti = 
#T2 = TO. Now let Ti — {ii, • • • , im}, ^2 — {j'l, • • • , jm} and associate to T the 
half-period characteristic fJT such that 



■•Ik 
fjip + niji^^Y^ f cj (4.2) 

7^ '^ 



where uj is considered as vector composed of first kind differentials. It can be easily 
shown that the the ambiguity in dividing of T into subsets and enumerating these 
subsets changes tJt by full period which can be ignored. It is convenient to measure 
the characteristics relative to Riemann characteristic S which is defined as 

flu, if 5 = l(TOod2) 
V{u\i}, a g = 0{mod2) 

where U = {1,3, 5,- • •}. The characteristic rjx corresponding to the subset T such 
that #T = g+ l{mod2) is defined by 

VT = f]ToU 

where T o U = {T U U)\{T n U). 

Our nearest goal is to prove the following 

Proposition. Consider the solution of KZ on level zero which corresponds to 
7ir'' i7ff = 0,1, ■ ■ ■ ,ag. The components of this solution are denoted by 
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/(Ai, • • • , A2g+2)'^^' !«2g+2^ ^q evcry multiindex ei, • • • , e2„ a partition of B cor- 
responds such that B = TUT', i e T : ei = +,i € T' : ei = -. Denote the 
elements of T by ip, p = 1, • • • , g + 1. The set T is associated with 9- characteristic 
rjT- The foUowing relation holds: 



/(Ai,--- ,A2s+2r^--'^^«+^ =C-30M(O)Met||9p9,log0M(O)|| 



ffXff 



(4.3) 



where C is the constant: C = rii<7(Ai ^ Aj)4A, A is the following important in 
future determinant: 



A = det|| / a. 

J ai 



jWgxg 



The notations are used: 

a,F(...,o,---) 



di 



-F{- 



To prove the proposition we have to explain certain facts concerning ^-functions 
on HES. 

Tomae formulae says that 6'[?7t](0) = if #r 7^ g + 1, and if #r == g + 1 then 



o[vTKof= n (^«-^.)^ n (^^-^.)^^ 

i<],i,jeT i<],ijeB\T 



(4.4) 



The following remarkable relation between the differentials u'^ the derivatives of 
6'-function at z = (6'-constants) takes place. If x, y, rj are such that 0[ri]{J^ uj) — 
then [20]: 

uj'{x,y) = -Y,^^^Jlog0[rlmuJ,{x)uJJ{y) (4.5) 



To proceed further we need explicit formulae for certain differentials. The nor- 
malized first kind differential LUi is given by g x g determinant: 



' ai(x), f CTl • • • f (Tl, f 



CTl, 



uj^{x) = (-1)*A-Met 



.^ai^)^ L,^9 ■■■ L_,^9^ /a.+i'^l 



9' 



L ^^ 



la, '^9 



Let us consider now the differential ll!'^{x, y) as function of x specifying y to one of 
the branching points: y = Xi. We take ^ = ■\/(r — Ai) as local parameter in vicinity 
of Ai. The differential u;^(x, Ai) is second kind differentil with zero a-periods and 
fixed singularity: uj'^{x, Ai) ~ ^~'^,x ~ Ai. It can be presented as {g + 1) x {g + 1) 
determinant: 



f'^ii^)^ la,'^^ ••• 



uj\x,X,)^A-\P'{X,))^det 



L ^1 \ 
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where Piir) = 



P'{Xi) = Yiji^i ^ ^j)- From this formula one 



derive the foUowing representation for uj'^{x,Xi) when x — Xj in terms oi g x g 
determinant: 



^(A.,A,)^A-(|;|M)^dct 



Jai ^9^3 



la, '^2,3 \ 

Jag ^9^3 
Lg P^,3 J 



(4.6) 



where aij = ai — Xjai-i, pij = pi + {Xj — Xi) ^u\. Combining (4.6) and (4.5) 

one arrives at 



Y.d^d^d\^TmAlAlXXX) = A-ip'(A.)det 



fla^'^^,3 

Jai ^9^3 
\ Ia^ P^J 



Lg '^2,3 \ 

Jag ^9^3 
lag P^,3 J 



(4.7) 

where i,j € T, the matrix A connects normalized first kind differentials with trivial 
ones: 



Lu == Aa, A 



ai 



(4.8) 



The reason for the equation (4.7) to hold is that adding J^' uj to r\T we get odd 
characteristics corresponding to T\{z, j}, and ^[77T\{i.j}](0) = 0. 

Let us return to the solution of KZ equations corresponding to 71 , • • • ,73 = 
oi, • • • , ag. It is given by (3.3): 



/(Ai,---,A23+2)^--'^^«+^-n(^.-A,)^ 

i<3 

n (A,-A,)-idet|| / Q{T\T\T')dT\ 



9^9 



(4.9) 



i£T,jeT' 



where T, T' are associated to ei, • • • , Cg in a usual way. Let us take g elements from 
T (say ip, p = 1, ■ ■ • , g) and multiply the determinant from (4.9) by Vandcrmond 
composed of corresponding A's the result being 



detllK 



f ^ 

"'"p m=l 



The following two differentials are equivalent (differ by a total derivative): 

,m^ig ^ Pq — XX \^iq^^3>~f^ ^ 

m=l j&B\T '^ ~ 

Hence the determinant in (4.9) can be replaced by 

I 






n (A* - A. 

i<i',i,i'£T\ig^i 



-det(C), Cj 



pq 



Pq 



(4.10) 
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Let us find the matrix X whicli satisfies tlic equation: 

AX^C (4.11) 

with A and C given by (4. 8), (4. 10). Kramer's rule tells that 



Xp,q = (-l)PA^Met 



la, '^P-1 
la, ^P+1 

la, '^9 



la, <^P-1 
la, '^P+l 

la, '^9 

Ln J 



Let us now consider the determinant of the matrix X . Multiplying it by Vander- 
mond composed of Ai^. , k — 1,- ■ ■ ,g one gets 



det(X) = 



1 



n iK-\') 



dct(X) 



where 



/L^2:F 



Xp,q = A"Mct 



L ^2,p \ 



Jai '^9,P Ja, '^9,P 

V /ai ^« Ia,^'1 J 



(4.12) 



It is easy to show that /i^ in (4.12) can be replaced by P'{Xi ) Pi ^i ■ Combining 
that with (4.7) and taking into account that dct(j4) = A^^ one gets: 

det(X)= Y[ [K - Xv)det\\dpdq\og6[T^T]m\gycg 

i<i' A,i'€T\ig^i 

Calculating determinants of RHS and LHS of formulae (4.11) and having in mind 
(4.8) we obtain: 

detll / Oil = Adet||ap9,log0M(O)|| (4.13) 

J ai 

Together with Tomae formulae it provides the representation for the solutions of 
KZ on level 0: 



/(A 



1; • • • 



.X2g + 2y^' 



.«2g + 2 



C-^0M(O)Met||ap9,log0[r?T](O)|| 



gxg 



where T is related to ei, • • • , £2^+2 as explained above, C = J^^^ (A^ — Aj)* A. That 
proves the proposition above. 

The formula (4.3) has very beautiful meaning. It relates the solution of KZ equa- 
tions which arc differential equation with respect to moduli of HES (the positions 
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of branching points) in terms of derivatives on the Jacobian. It would be nice to 
prove directly that (4.3) satisfies KZ using the heat equation for the ^-function. 

Now suppose that we took a solution corresponding to other choice of half basis: 
7i, • • ■ ilg — '^'ii' ' ' io!g- Denote this solution by /^^/(Ai, • • • , X2g+2)- Clearly, all the 
reasonings above are applicable to Ja' the difference being that the final formula 
has to contain 6*- functions defined with respect to the half-basis A' . In order to 
rewrite the answer in terms of canonical 6'-functions (those corresponding to A) we 
have to use the formula (4.1). After simple computations one gets: 

/a'(Ai,---,A2,+2)^--'^^«+^ = 

det(M)C-30[77T](O)4detPp9,log0[r;T](O) + ^,.^l°g(detM) ,, 



where 



i)-{'rMl.Y it 



d^ij 



^ ; eri.2. 

a 



The matrix M is defined as M — {cfl + d) . 

In Section 6 we shall need another representation for the solutions of KZ in terms 
of 0- functions. The author was unable to prove this representation completely, so 
we shall formulate it as conjecture and present reasonings in favour of it. 

Conjecture 1. For every half-basis A' a polynomial Q^i in di of total degree 
2g exists such that 

fA'iXl,--- ,X2g+2r''- ''''+' =C-''e[TjT](OfQA'0[VT]{Of (4.14) 

where the connection between ei,- • • ,£29+2 and T ia as usual, the coefficients of 
Qa' might be complicated, but they do not depend on T. 

Let us explain why we assume the representation to exist. To this end we have 
to understand what kind of constants depending on T is a candidate for being 
presentable in the form 9[riT]{0) Q9[riT]{0) for some polynomial in di and T corre- 
sponding to even non-singular characteristics (#r = g + 1). Riemann ^-functions 
on HES satisfy many relations which follow from the combination of Riemann re- 
lations with the peculiar properties of HES. In particular, the following relations 
hold for 9[ri]{z) (rj is even, nonsingular, i.e. such that 6'[r7](0) 7^ 0)[19]: 

E {-lf''^''°^^e[rjs]{zf9[7^s]{0f = (4.15) 

ScB,#S=g+lAeS 

where T is arbitrary subset of B satisfying the requirements :#T =17-1- l{mod 2), 
=ffT =^ g+l, 1 € T. We put the requirement 1 € 5 in order to avoid the summation 
over B\S. So, if the set of constants is presentable in the form Q9[rjs]{0)^ then 
they should satisfy this system of relations being placed instead of 9[ris]{z)'^. The 
opposite should be also true: if a set of constants enumerated by S satisfies the 
system (4.15) then they can be presented in the form Q9[ris]{0)'^ ■ In elliptic case 
(g — 1) there is only one set T satisfying the requirements above: T = B. For this 
set the relation (4.15) turns into 

eo.o(O)'0o,o(^)' = eism^9^,^{zf + 0o,i(O)'^o.i(^)' 
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in usual notation for the elliptic ^-functions with characteristics [19]: 

^1,0(2) == ^['7{i,2}](2), ^0,0(2;) = 0[V{i,3}]{z), ^0,1(2) = 0[?7{i^4}](z). 

Let us prove that the solutions of KZ for arbitrary choice of A' satisfy the 
relations (4.15). First, let us check that 6'[775](0) given by Tomae formulae satisfy 
the requirements. Substitute 9[ris]{0) given by (4.4) into (4.15) instead of both 
^[^s](0) ^^^ ^[Vs]{^) s-nd divide the relation by A^ ni<7(^« ^ ^j) ^^'^ result being 

E (-ir^^°^^ n a -A^ -Q (4.16) 

These identities are proven by consideration of residues at Aj = Aj, Vi, j. 
Now let us present the solution of KZ (4.9) as multidimensional integral: 



fA'{\l,- ■ ■ ,\2g+2r'-^''^+' =l[iX^ - Xj)' 



X n {X,-\,)-' f dT,^l=--- j dT,^l=dct\\Q,{T,\S\S')\\,^, 

where Qi are given by (3.4), S is related to ei, • • • , £2^+2 as usual. Let us denote 
det||(5i(Tj)||gxg by Xg{Ti,- ■ ■ ,Tg\S\S') indicating the dependence on genus g and 
subsets S,S' (recall that S" — B\S). We shall prove that the solutions of KZ 
satisiy (4.15) being substituted as 6'[775](0) d[ris]{z) if we prove that Hies jes'l-^i" 
\j)~^Xg{Ti, ■ ■ ■ ,Tg\S\S') satisfy them. So, we have to prove that 

^ (_l)Sn(aoT) -Q (A,-A,)-%(ri,---,r,|5|5') = (4.17) 

scB,#s=g+iAes les.jes' 

In fact the proof of the identity does not differ much from the proof of (4.16). 
Again it is sufficient to check the cancellation of the residues at A^ = A^ because 
^g{'''ii ■ ■ ■ T^glSlS') is of degree g with respect to any A. Suppose I Cz S, m € S' 
then the polynomial Xg{Ti,--- ,Tg\S\S') satisfies the following recurrent relation 
(classical version of the relations from [11]): 

X,(ri,...,r,|^|^')L,=,„-n(r.-AOE(-l)'{^ IT i^, ~ X,)} 

k=l fe=l '' peB\{l,m} 

X Xg_l(Tl, • • • ,Tk-l,Tk+l, ■ ■ ■ ,Tg\S\l\S'\m) 

The point is that the coefficients in the relation are independent of 5*. That provides 
the possibility of inductive proof of (4.17). 

It would be nice to have explicit formulae for the polynomials Q in (4.14). The 
author has not succeeded to get them. Presumably the formulae can be found from 
(4.3) using Fay identities [20]. 

To finish this section let us present explicit formulae for the elliptic case. The 
canonic way for construction of homology basis is as follows: the cycle a surrounds 
the cut between Ai, A2, the cycle b starts from upper bank of this cut, reaches 
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the upper bank of the cut A3, A4 by one sheet, then moves to another sheet and 
starting with the lower bank of A3, A4 returns to the lower bank Ai, A2. There are 
two independent choices of a-cycle: a — a or a' — b which are connected by the 
following clement from Sp{2, Z): 

' 1 
-1 

So, there are two independent solutions to KZ: 

/a(Ai,--- ,A4r'-''^ = c-30M(o)*a2iog0M(o), 

(4.18) 
A2 

where K is full elliptic integrah/iT = i / .'^"^ ,the correspondence between ei, • • • , e4, 

" Ai V ^f'^) 

T and conventional notations for elliptic ^-functions are as follows: 



{++—},{- 


--++}- 


-{1,2}- 


-{1,0}, 


{+--+},{- 


-++-}- 


-{1,4}- 


-{0,1}, 


{+-+-},{- 


-+-+}- 


-{1,3}- 


-{0,0} 



(4.19) 



The solution can be also rewritten as (4.14) with 
since dO[rj\{Q) ~ for even -q. 
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5. Pragmatic view of finite gap integration. 

Consider an integrable equation with infinitely many times ti,t2,- ■ ■ ■ To any 
time ti ttie M-operator dJli is attaclied. Tlie M-operators satisfy zero curvature 
condition: 

We consider s^(2)-invariant case (one can think of nonlinear Shrddinger model for 
example), so 971^ is traceless 2x2 matrix depending on spectral parameter A such 
that 

i 
tTm^{\) = 0, m^{X) = Y^ X'^rUi^k, 
fe=0 

also we normalize it by the requirement ruii — a^ . 

The coefficients of the M-operators are dynamical variables. Finite-gap integra- 
tion [16-18] deals with the situation of stationary solutions which means that there 
is a time <„ on which the dynamical variables do not depend: 

J-OT, = 0, Vi 

That means that the M-operator 9Jl„ satisfies the equation: 

hence the determinant of OJln is an integral of motion with respect to all the times. 
Certainly this determinant can be presented as follows: 

In 



detaJl„(A) = [](A-Afe) 



fc=i 

That suggests that the problem in question is closely related to the problem of 
parametrization of all matrices M(A) (we omit the index n) which satisfy the 
following requirements: M(A) is traceless 2x2 matrix depending on A as poly- 
nomial of degree n with fixed senior coefficient (equal tr"^) and given determinant 
(detM(A) = nfc=i(^ ~ ^fe))- This problem was considered in the last century (Ja- 
cobi, Riemann), the solution can be found for example in Mumford's book [19]. Let 
us describe it in the terms appropriate for our further goals. 

The matrix in question is degenerate at the points A = Aj, also it is traceless, 
hence M(Aj) is a Jordan cell: 

for some vector ipj. Let us parametrize M(A) by the set ^j, j = 1, • • • , 2n (these 
vectors are not independent as we shall sec soon). Construct the interpolation 
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Definitely, M' is degenerated at the given points with given values, but its degree 
in A equals 2n — 1 instead of n required. So, we have to kill n — 1 senior coefficients 
of M' , also we have to take into account that A" is supposed to enter with given 
coefficient a^ . After some simple manipulations these requirements lead to the 
following system of quadratic relations for the components a^, hi of the vectors f/^i: 



y 1 

^^ rr' (\ 



keT Ujeri^k - A, 



-4 = 0, 



y I 

^-^ rr' f\ 



k£T '■'■] 



Yljeri'^k - Aj 



-bl = 0, 



fcgTlljGT(Afe-Aj) 

for any subset T of n+1 elements of the set {1, 2, • • • , 2n}. This system of equations 
leaves only n — 1 independent parameters (which could be, certainly, calculated from 
the very beginning). The important point is that due to Ricmann identities and 
Tomae formulae the relations allow a parametrization in Riemann 9 functions on 
the hypcrcUiptic surface r^ — P(A) = 11/0=1 (^ ^ ^k)- The solution to the system 
(5.2) looks as follows: 



9[S]{z)9[rj,]{ry ■' e[5]{z)e[rj,]{r) 



^ C ^'fZr/^ ^ h = C JcfClr u\ (5-3) 



where rij is the theta-characteristic corresponding to the branching point \j : r]j — 
f]\^j with f]T defined by (4.2). The variable z = (zi, • • • , Zg) is the parameter on the 
Jacobian (genus g = n — 1), 

/■oo"*" 

r^ 9 tj, 

00+ is one of two infinities on the surface, u is the vector composed of first kind 
differentials, S is Riemann constant, finally, C is the same as in (4.3). 

We would like to make two comments on the above formulae in order to clarify 
them from two points of view. First point is the connection with more familiar in 
the context of finite-gap integration object, namely with Baker-Akhiezcr function. 
The BA function is an eigenvector of M(A): 

M(A)V'(A) = m(A)V'(A) 

where m is corresponding eigenvalue. Certainly, for generic A there are two solution 
to this equation which means that the BA function is defined on the surface S. But 
for A = Aj the matrix M{X) is degenerate, it has only one eigenvector corresponding 
to zero eigenvalue. This eigenvector is exactly our ipi. Thus, 

ipi = tpiXi) 

The BA-function tp{x) {x = A* is a point on the surface) is written in terms of 9- 
functions through 9(J lu + ■ ■ ■), [16-18] when x coincides with one of the branching 
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points this expression turns into 0-function with corresponding half-period charac- 
teristics because the integral of lo taken between two branching points on hyperel- 
liptic surface is a half-period. 

The second point is the direct connection with well known addition theorems for 
0- functions in the elliptic case (n = 2, there are four branching points). In that 
case the following identification can be done with usual ^-functions [19]: 

0[r;i](z) = 0o,oW, 0[r;2](z) = eo,i(z), 

Also the following simple variant of Tomae formulae holds: 

00,0(0)2 ^ [(^^ _ y^-^^y^ _ x^y^h^^ 0p^^(o)2 = [(Ai _ A4)(A2 - \^t-K 
^i^o(O)' = [(Ai - A2)(A3 - U)\^K, 0ia(O)2 = 



where K is the full elliptic integral K — -k J dr/ ^yPij), a-cycle surrounds the cut 
between the points Ai and A2. The ratios of 0- functions at r are easy to calculate: 

0fe](r)2 = Const P'iX.rK P'{\,) ^ \{{^j - Afe) 

Taking all that into account one realizes that for the paramctrization (5.3) consid- 
ered the equations (5.2) turn, for example for T = {1,2,4} into 

Oofii^fOos^r + zf - 0o,i(O)'0o,i(r + zf - 0i,o(O)20i,o(r + zf = 0, 
Bofii^feoAr - zf - eo,i(O)'0o,i(r - zf - 0i,o(O)'0i,o(r - zf = 0, 
dofii^fe^^r + z)9oAr - z) - 0oA{Ofeo,i{r + z)0o,i(r - z)- 
OiAOfOi,o{r + z)9iAr - z) = di^izfOiArf 

The equations (5.2) for other subsets T produce other known identities for 9- 
functions. In the case of generic n the situation is similar: in paramctrization 
(5.3) the equations (5.2) turn into certain special cases of the Riemann identities 
(Frobenius formulae [19]). 

Returning to the integrable models we conclude the following: the formulae (5.3) 
provide the paramctrization of Af (A) on the torus ( Jacobian), that is why the times 
U should be related to Zj as 

9 

the constants c^ refer to the particular integrable equation which is not essential 
for us. 

Let us consider now the same situation from different point of view. What will 
follow is an extraction from Adler-Reyman-Semenov-Tian-Shansky approach[21,22] . 
Consider the loop algebra g — sl(2) with generators J^, a ~ 1, 2, 3; m is an integer. 
The generators satisfy the relations: 
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The algebra contains two subalgebras: 0_ generated by J^ with to < and g+ 
generated by J^ with to, > 0. Evidently, the algebra allows finite-dimensional 
representation py. 

By G = SL{2) we denote a groop whose Lie algebra coincides with g. This group 
contains two subgroups: 

G+ = cxp(0+), G- = exp(0_) 

Consider now the matrix M(A) of the same type as above (traceless, polynomial 
in A of degree n with fixed senior coefficient) and introduce the action of G on M(A) 
as follows: 

g{M{X)) =. px{g)M{X)p^{g)-\ for g e G (5.4) 

The subgroups G^, G- acting on M{X) generate the orbits 0+, O- of the form: 

oc 

ff+(M(A)) = ^TO;,A^ 

fe=0 
n 

g-{M{\))^ Y, <^\ < = f^3 

fc— — oo 

The algebraic interpretation of the integrable models is explained as the problem of 
description of the intersection J of the orbits Oj^^O-. Clearly this intersection is 
composed of the matrix of the same type as M[\) itself, so the connection with the 
matters discussed above is manifest. How to describe the intersection in question 
in algebraic way? If M' E J then there exist g+ , g_ such that 

g+{M{\)) = px(g+)M{\)px{g+)-^ = 
= g-{M{\)) = px{g^)M(\)px(g^)-^ (5.5) 

It is evident from these equations that the matrix .g(A) = p\{g^)^^ p\{g+) commutes 
with M(A). The only possibility for the matrix to commute with M(X) which is 
supposed to be not degenerate for generic A is to be a function of M(A). Fortunately, 
additional dependence on A makes the situation nontrivial. We can consider g{X) 
of the form: 

.g(A)=exp{^t,Af,(A)} (5.6) 

where Mi (A) = A-*M(A). 

Now let us inverse the reasonings which means to start with 5(A) in the form 
(5.6), and to try to construct J. Evidently to this end we need to solve the following 
Riemann problem: present 

g{\) = exp{^t,M,(A)} = px{g-)-^ px{g+) (5.7) 

where the Loran series for p\{gj^)^ {px{g-)) contain only positive (negative) powers. 
It is clear that only those Mi are essential for which < i < n, others can be directly 
moved to either g^ or g_, and do not contribute to (5.5). So, we deal with the 
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dependence on n — 1 times which can be, actually, identified with the parameters 
on the Jacobian considered above. 

Suppose the problem (5.7) is solved. Then we can introduce "times" dependent 
matrices M(A, ii, • • • , ig) as being dressed according to (5.5) by g+ or g^ from (5.7). 
ft easy to show that the projections of A~*Af (A, ii, • • • , ig) onto g+ denoted by 2Hi 
satisfy the equations: 

which coincide with (5.f). 

We would like to emphasize the importance of the vectors ipi from the algebraic 
point of view. The BA function itself does not transform in a reasonable way under 
the transformations (5.4), but its values at the branching points {il^i) do transform 
under finite-dimensional representation of G when M transforms under (5.4): 

ipi -> G{Xi)'ipi 

This nice property of ipi will be important in what follows. 

To finish this section let us write down explicitly the equations in terms of mf: 

^K= E rX-«+feS (5-8) 

q—max[0,k-\-j — n) 
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6.Tau-function and KZ equations. 

The most mysterious object in the theory of classical integrable equations is that 
of T- function [23-26]. Hirota observed that the integrable classical equations can 
be rewritten as follows. Let us consider certain function t, depending upon the 
arguments t — {ti, t2, ■ ■ ■} and introduce the notations: 

A,---A.r.r(t)^^...^r(i + x)r(i-x)|_„ (6.1) 

Then according to Hirota the integrable equations can be written as 

P{D)t-t = (6.2) 

for certain "polynomial" in D. We shall be interested mostly not in t- function 
itself, but in the linear space of functions 23 which is defined as a space spanned by 
the functions [Di-^ ■ ■ ■ Dif,T ■ T{t)]/T{tY . Hirota equations (6.1) imply certain linear 
dependence in the space. 

The main achievement of Kyoto school in that direction is the formula expressing 
T- function in terms of vacuum expectation in the space of highest weight representa- 
tion for the algebra sl{2) on level 1 [24]. Consider the Cartan subalgebra generated 
by /fe = J|, fc > 0. Then the formula takes place 



T(f) = (r!|.gexp(^ife/fc)|r!) 



where [fi) is the highest vector of the representation on level 1, {Vt\ is the dual 
vector, g is the element of central extended group SL(2) which specifies a particular 
solution of the equation. We want to relate the KZ equations on level zero to r- 
function. The formula (6.2) is both good and bad for our goals. Good thing is that 
(6.2) links the integrable equations to the highest weight representations which are 
involved in KZ. Bad thing is that the formula (6.2) on the one hand and our KZ on 
another deal with different central extensions of sl(2). We shall take the good thing 
as a hint of possible connection between fc = KZ and r-function, and forget about 
the bad thing. It should be said also that for the case we are interested in (that of 
finite-gap integration) the description of g is rather complicated and indirect. 

What is T-function in the finite-gap case? The answer is well known: it is 
essentially the 6'-function 6'[(5](z), ti = '^{.c'^Zk, where c are some coefficients. 
How is it connected with BA-function? The usual answer is that it coincides with 
the value of BA-function at certain point (oo). This answer does not satisfy us 
because we do not like to use BA-function at generic point for its unclear algebraic 
properties, we want to deal with the values of BA-function at the branching points 
(ipj) only. 

The space 2J is finite-dimensional in finite-gap case. Let us consider the example 
g — 1. All the times U are proportional to one variable on the Jacobian z. The 
T-function is di^i{z). The space 03 is generated by 1 and ^^l''^ ~ ■^\og9i^i{z). 
The Hirota equation reads as 

D'^T ■ T + CiDV • T + C2t'^ = 
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which turns into well-known differential equation for 0-function (equivalent to the 
equation for the Weierestrass *P- function) : 

^A ^2 f)2 

^log0i,i(z) + 6(— log0i,i(z))2 + ci—logOi^iiz) + C2 = 

for some Ci,C2. 

Our further strategy will be in putting together the space 2J, the vectors tpi a-nd 
the solutions of fc = KZ equations the idea being the following: Vi transforms 
under finite-dimensional representation of sl(2), r is connected with highest weight 
representation, hence they should be put together via KZ equations. To start let 
us consider the following object: 

it should be said that v{z) is not a function but rather a set of functions because 
different choices of blocks of vertex operators (intermediate Verma moduli) are pos- 
sible. This is the same as different choices of cycles A' for /^/(Ai, • • • , X2nY^'"' ''^^"'■ 

{0\V''{Xl) ■■■V''-{X2n)\0)so,ne block = ./.omeA' (Al , • • • , Aa™)'^ ■- ■'^" 

So, if necessary we shall indicate a particular component of v as va' ■ The first 
important property of the functions v{z) is that their dependence on z is governed 
by "free" dynamics in Verma module. 

Let us take the matrix M(A) as described in the previous section: 

n 2n 

M{X) = J2 A'"^>"- < = Sa.^3, det(M(A)) = [](A - A^) 
Associate to the matrix M{X) the clement of g 



4=0 

such that M(X) = p\{M). Similarly, the algebra elements Mj can be introduced: 

n 

M, = J2 <Jt-j, Px{M,) = M,{X) 

i=0 

Clearly, Mj commute among themselves: 

[M„M,]=0 

Let us consider M{X) as a starting point for the procedure of the previous section, 
i.e. M{X) is the stationary M-operator corresponding to the moment tt = 0. To 
the matrix M(A) we attach the set of vectors '^|}^^ i = 1, • • • , 2n as it is explained 
above. The following important statement holds. 
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Proposition The z-dependcncc of v(z) is governed by "free" dynamics in the 
Vcrma module: 



viz) = ^i^,,{z) ■ ■■i^2n,.,AzmV"'{^l) ■ ■ ■ V''"iX2nm = 

n-1 

V5i,.x(0) • ■■i>2n.e,A0){0\V''{Xi) ■ ■ ■ F^^"(A2„)exp(^ M,U)\0) (6.3) 

i=l 



where z and t are connected by some hnear transformation: U = c^^j. 

In what follows we shall allow ourselves inaccuracy using both notations: tpiiz) 
and ipiit) assuming that they coincide when z and t are properly related. As it was 
explained in the previous section the functions ipi satisfy the following equation: 

d - 

--^^it)^Mt)MK,t) 
atk 

where the t-dependence of 971 is due to (5.8). This equation allows to express all the 
higher derivatives of ■0 at the moment t = in terms of derivatives of $H. Let us 
prove the equations (6.3) comparing all the derivatives at t = of LHS and RHS. 
With first derivative it is trivial, in LHS one has: 

^(ViM,(i) • • •^2„,..„(i)(0|F'HAi) • • • V^^^"(A2„)|0))|i.o - 

2n 

= ^ 1^1,,, (0) • • • vip-i,.,_, (0)v5p,,. (o)^p+i,.,+, (0) • • • 7^2„,..„ (0) 
p=i 

In RHS one has: 

o n— 1 

■^i>l,.^{0) ■ ■■^2n,e,A0){0\V''{\i) ■ ■ ■ F^^ (A2„)exp(^ M,U)\0)\t=O = 

= ^l,.i(0) • ■■^2n,e,A0){0\V'H^l) ■ ■ ■ F'^"(A2„)Mi,+ |0) - 
2n 

= J2 ^1.^1 (0) • • • V^p-1,.,-1 mp,e', (0)V^p+l,e,+ i (0) • • • V^2„,e.„ (0) 
p=l 

x<^ (Ap)(0|F^HAi)---T/^^"(A2„)|0) 



Evaluating the RHS we used the notation Mi^+ for the projection of Mi onto 0+; 
recall that J^|0) =0, m < 0. We moved Mi^+ to the left using the properties of 
the vertex operators, Mi,+ annihilates the left vacuum, also Vyii{\p) ~ p\ (Mi,+). 
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Less trivial computation is that for the second derivatives. In LHS one gets 

= J2 ^i,.i (0) • • • v^p-i,.,-, mp,e'^ mp+i.e,+, (0) • • • 

p<q 
X V'q-l^e,-! {0)i^q,e'^ (0)V'g+l,e,+ i (0)V'2n,e2„ (0) 

X {Ti^^^,,{Xp)m%^{x,) + Ti^^^,,{x,)m%^{x,)}{o\v^^ (x,) ■ ■ ■ v^^-{x,„)\o)+ 

2n 

+ Y, ^1,,, (0) • • • vip-i,.,_, (0)v5p,,, (o)Vip+i,.,+, (0) • • • vi2„,..„ (0) 
p=i 

X {OT:i;„(Ap)9}l<^(Ap) + 9,9?i;>,^(Ap)}(0|y'KAi) • • • T/^^"(A2„)|0) 

To transform the RHS one has to deal with MiMj\Q) which can be evaluated due 
to the equations (5.8): 

M,Mj\Q) = (M,,+ + M,,_)M, |0) = 
= Mi,+M,- +|0) + a,Mj,+|0) (6.4) 

Now moving Mi^+Mj_+ and diMj^+ to the left one gets (6.4). The consideration of 
general case is quite similar to the case of the second derivative which is the most 
demonstrative one. 

Thus, we have shown that the functions v{z) possess the nice property (6.3) 
which makes them similar to T-function. 

It should be mentioned that the starting M{X) can be taken in especially simple 
way. For example we can divide the set of Xi into two subsets {B = SiJ S', #5* = 
#5" = n) and define M as follows: 

/n(A-Ao+ n(A-A,), n(A-A.)- n(A-A. 

M(\\ — I «e-S «es' ies ies' 

^ '~ n(A-A.)-n(A-A,), -n(A-Ao- n(A-A. 

ViGS' les ies les' 

In that case the vectors tpi are also simple. 

We have to connect the functions v{z) with some familiar objects. There should 
be a clever way to get the result which will be announced soon. But we do not know 
that way which makes us to proceed to calculations. Let us consider the elliptic 
case. Recall that the covectors ^'i are given in the elliptic case by: 

J , N /^o.i(?' + z) 6*0, i(r - z). 
tp2(z) ~ (- ■ 



Mz)^C 



lAr + z) Oi^r-z), 



0i,i(z) ' 0ia(z) ^ 

T , -, fOi,o{r + z) 6'i,o(r-z). 

Sia[z) oi^i{z) 
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In what follows we shall perforin calculations up to constants which will be outlined 
by using ~. There are two solutions to KZ equations in elliptic case which are given 
by are given by 

the correspondence between ei,--- ,£4, T and conventional notations for elliptic 
0- functions is given by (4.19). This formulae suggest to generalize the expression 
for v(z) substituting in it 9[riT]{0)'^d[riT]{x)'^ instead of the solutions of KZ. Doing 
that one gets the following: 

v{z,x) ~ 6*1^1(2)^ 

+ Oo.oir - z)eo.i{r - z)6ii,i(r + z)6li,o(r + z)}+ 

+ OoAOfOo,i{xf{eo,oir + z)9oAir - z)0ia(r - z)0i,o(r + z)+ 

+ 9oAr - z)9oAr + z)6'ia(r + z)eifi{r - z)}+ 

+ OoAOf^o,oixf{^oA^ + ^)^o,i(^ - ^Whiir + ^)^i,o(^ - z)+ 

+ Oo.oir - z)9o.i{r + 2)6*1, i(^ - z)ei,o{r + z)}] 

This expression can be simplified via Ricmann identities the result being quite 
simple: 

/ X OiAz + x)di^i{z - x) 

V{Z,X) ~ —2 

which leads to the following nice formulae for Va{z),Va'{z): 

D'^T ■ T D^T ■ T 
Va{z) ^ 5 {z), Va'{z)~K ^ (2) + 1 

where d\^\ is denoted by t, D is Hirota derivative. 

Thus by convoluting -0^ and the solutions of KZ in the elliptic case we got exactly 
the basis of the space 03. Combining that with the proposition above we also realize 
that the dynamics in the space 5J is governed by the free dynamics in the Verma 
module. We would suppose that this nice connection holds in generic case as well. 
In order to prove that we have to use the Conjecture 1 from Section 4 together with 
the following 

Conjecture 2. The following identity (up to neglectable constants) for the 
^-functions on HES holds: 

2g+2 
T(lB,#T=g+l j=l 

c::i e[6]{zf^e[6]{z + x)9[S]iz - x) 

where T and {ci} are connected in usual way. 

The author strongly believes in this conjecture, still he has been unable to prove 
it. 
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If the conjectures 1,2 are both true in generic case then the following relation 
should take place: 

21 3 9^1-^ = ^ . . . v;2„,,,„/a,(Ai, • • • , X2nY''- ■''- (6.5) 

where r(z) = 6'[(5](z). We would suppose also that there are enough 9- constants 
hidden in /a', i.e. that all the independent Hirota derivatives of r can be obtained 
by linear combinations of (6.5) for different A' . 

The formula (6.5) establishes the relation between ^i which transform under 
finite-dimensional representations of sl{2) and the space of Hirota derivatives of 
T-function (the space of classical fields) via the solutions of KZ equations on level 
which we can call classical form factors. This is exactly the formula we were looking 
for, its similarity with quantum formula (1.1) is manifest. 

To finish this section we would like to present one more interesting formula: 

{0\V'' (Ai) • • • y^- (A2„)M'"exp(^ t^M - 

= V-i^O) • • • V'2^"(0){V5i.<(0 • • •^2„,4„(i)(0|H(Ai) • • • y^-(A2„)|0)|t.o} = 

= vr(o)-->2?r(o)«(t) 

This formula is due to the fact that 

PA.(M)2 = (V;,(5§V;,)2 = 
Notice that as usual M and Ai,--- , A2„ must be connected: det(pA (.^-^)) = 

n(A-A.). 
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7. Conclusions . 

Let us return to the very beginning of the paper. We were wondering whether 
a formula similar to (1.1) exists in classics. We have answered this question in the 
last section of the paper. Really, the formula (6.5) is very much similar to (1.1). 
It connects classical local fields (Hirota derivatives of r-function) with the tensor 
product of finite-dimensional representations of tha affine algebra via the classical 
form factors (the solutions of KZ on level 0). So, the similarity of the formulae 
(1.1) and (6.5) is described by: 

4>,, ^P, ^ Z{(3i), Z*{l3i) 
53^ 21 
f ^ F 

The quantum formula needs also the summation over all particle states which can 
be interpreted as summation over all HES, the number of particles being the re- 
lated to genus of surface, the rapidity being considered as the positions of branching 
points, i.e. the moduli of the surfaces. Notice that in the classical limit the dif- 
ference between the finite- dimensional representation and conjugated one (which 
were responsible for particles and antiparticles in the quantum case) essentially 
disappears. 

So, we would suppose the following procedure of quantization to be possible. 
We start with the family of finite-gap solutions of integrable model with sl{2) 
symmetry. Different solutions are parametrized by the intersections of the orbits of 
0-1- as it is explained in the Section 5. These intersections are parametrized by the 
matrix M(A), but we prefer to parametrize them by the set {V'i} connected with 
M(A) as it is explained above. The vectors f/'i transform under finite-dimensional 
representation of the affine algebra. We can take 



oo .. 

^c/-0 / Vx,®---®Vx^^ (7.1) 



""^Ai<---<A2„ 

{V\ is the space of the representation p\) as a, completion of the manifold of classical 
finite-gap solutions. Not every vector from this space is good for the classical 
solutions but only those which can be presented in the form ■01 (g) • • • (g) V'2n with 
the set {V'i} associated to a matrix M(A) as described above. For every particular 
element of this form the classical fields are given by convolution with classical form 
factor. In fact the space (7.1) coincides with the space of states of quantum model 
Hp (1.3). Also the quantum form factor is a quantization of the classical one 
which is essentially due to the quantization of the algebra underlying the theory 
{sl{2) —^ D{Y)). So, we propose that the following way for quantization should be 
possible: we extend the set of classical solutions to the space Hd and identify it with 
the space of states of quantum model, the rest of quantization is in quantizing of 
symmetry algebra. It should be strongly emphasized, however, that the space-time 
of the quantum model has nothing to do with the classical "times" ti, it appears 
as a result of quantization as it is explained in Section 3. 

Let us finish with several remarks. The finite-gap solution constitute in the 
classical theory nice but small subset of solutions. Generally, there are infinite- 
gap solutions in which the finite-gap ones are of measure zero. The infinite-gap 
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solutions are rather ugly ones, no reasonable theory is available for them. So, since 
for the quantization only finite-gap solutions are needed, we have a nice example of 
usual phenomena: quantization takes everything good from the classical theory and 
forgets about bad things typical for it. There is a puzzling connection between the 
matters discussed in the present paper and those considered in [27] in the connection 
with the theory of strings (this remark is due to H.Oogury). The constructions of 
the present paper can be generalized to su{N) case for which the algebra sl{N) 
and the Riemann surfaces with the branching points of A^-th order are responsible, 
the form factors of corresponding su(A^)-chiral Gross-Ncveu are given in [11]. But 
probably this is not the best way for the generalization. We would better proceed 
to the consideration of arbitrary Riemann surfaces which are associated to KP- 
equation. This case should cover all the su(A^)-invariant Gross-Neveu models and, 
probably, will lead to something essentially new. 
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